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ABSTRACT 


Intensity oscillations of coronal bright points (BPs) have been studied for 
past several years. It has been known for a while that these BPs are closed 
magnetic loop like structures. However, initiation of such intensity oscillations is 
still an enigma. There have been many suggestions to explain these oscillations, 
but modeling of such BPs have not been explored so far. Using a multithreaded 
nanoflare heated loop model we study the behavior of such BPs in this work. 
We compute typi cal loop lengths of BPs using potential held line extrapolation 
of available data fjChandrashekhar et al.ll2013l) . and set this as the length of our 
simulated loops. We produce intensity like observables through forward modeling 
and analyze the intensity time series using wavelet analysis, as was done by previ¬ 
ous observers. The result reveals similar intensity oscillation periods reported in 
past observations. It is suggested these oscillations are actually shock wave prop¬ 
agations along the loop. We also show that if one considers different background 
subtractions, one can extract adiabatic standing modes from the intensity time 
series data as well, both from the observed and simulated data. 


Subject headings: methods: data analysis, methods: numerical. Sun: corona. 
Sun: hares. Sun: oscillations. Sun: UV radiation 


Introduction 


Coronal bright points (BPs) are ubiquitous fea tures of the q u iet-Su n and of the coronal- 
hole regions. These features were hrst sighte d bv Ivaiana et ahl ( 1973 1 while observing the 
Sun using rocket born X-ray telescope. Later, iGolub et al.l (1197411 determined the lifetime of 
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such X-ray BPs to be about eigh t hours using Sk ylab X-ray images. BPs are also obsery ed in 
EUV f Habbal &: Withbroe 1981 ). It was found ( Krieger et ah 1971 : Golub et ah 1974 ) that 
these BPs are actually small magnetic loop like structures associated with mag netic bipolar 


regio ns in the photosphere. The typical size of these BPs are about 10® Km^ flGolub et ah 

19741 ) ■ 


The most distinctive feature of these BPs is their temporal variability during their life 


time 

Jsheelev & Golub 

1979; 

Habbal & Withbroe 

1981; 

Kumar et ahlboil 

Ghandrashekhar et al. 

2013 

). It has been suggested that this temporal variability is either due to the interactions 


of the solar p-mo de oscillations with t 
netosonic wav es flBogdan et al.l l2003l: 


reconnections flDoyle et al 


2006 


le magnetic fields in the corona resulting into mag- 


Kuridze et ahl 120081) . or due to recurrent magnetic 


Tian et ahl 120081) . Our primary aim in this work is to in¬ 


vestigate the extent to which the latter possibility contributes to the intensity oscillations of 
the BPs. 


Since BPs are composed of small compact loops (IGolub et al.lll974lKwon et al.ll2010[) . it 
will be sufficient to numerically simulate a small loop system to model coronal BPs. It is well 
known that due to low plasma beta of the lower corona, plasma is conhned along the coronal 
magnetic held lines. This enforces the plasma to show one dimensional (ID) dynamics along 
the held lines, allowing us to model the coronal loop in ID. While observing through space 
or ground based telescopes, however, it is not possible to resolve every individual strand of 
an observed coronal loop. Rather, what we observe as a global loop is actually an ensemble 
of several ID stra nds together. It is universally accepted that swarms of tiny reconnections 
called nanohares fjParkeii 119881) . each releasing energy around 10^® ergs, are responsible for 
the high temperature (above 1 MK) of coronal loops. 

The purpose of this work is to see if such nanohare like heating events can give rise 
to density huctuations that are rehected in the high density chromosphere before bounc¬ 
ing b ack to the corona again, thereby emulating the observed oscillations in the coronal 


BPs flKumar et ahl l201ll: iGhandrashekhar et al. 120131). To tha t end, we employ the ID 
hydrodynamic coronal loop code developed in ISarkar fc Walsh! (120081) . To replicate the 
multi-stranded nature of the global coronal loop, we simulate 125 ID hydrodynamic strands 
individually. These simulated strands are heated by nanohare like heating events. Since 
BP loops are very tiny (~ 10 Mm, see below) compared to active region l oops, we have 
assum ed that loop expansion ( Mikic et al. 2013 ) or thermal non-equilibrium ( Lionello et al. 


20131 ) does not play a crucial role in the dynamics of the loops. From the simulated global 
loop, we produce observables like intensity, and hence deduce the oscillation periods to see 


if they agree with the real oscilla tions studied through various instruments flKnmar et ah 
2011 : Ghandrashekhar et al. 2013). 
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The r est of the paper is organized as follows. In Section |2] we briefly discuss the obser¬ 
vations of IChandrashekhar et ahl (120131) and from the observed BPs we hnd out the typical 
loop length which we take as input of our simulation. In Section [3] we describe the simula¬ 
tion details of this multithreaded coronal loop model. In Section H] we describe the forward 
modeling to produce loop intensity in the AIA bands, which we analyze in Section [5l Finally 
we conclude in Section |6] with a summary of this work and proposal for future observations. 


2. Observed intensity fluctuation and loop length determination 


In IChandrashekhar et al.l (120131 ). ten different BPs were observed using full-disk Level 1.0 
AIA and HMI images obtained over a period ranging from February 13 to February 15, 2011. 
It was found that each BP intensity varies with time. Wavelet analysis on every BP image 
was also performed to measu re the prominent periodicitie s. Further details of this work 
can be found in section 2.1 of IChandrashekhar et al.l (120131 ). In the present work, we want 
to compare BP intensity period icities computed from our simulations with those measured 
bv IChandrashekhar et al.l (120131) . 


As mentioned earlier, we plan to simulate a multistranded small coronal loop which can 
represent a typically ob served BP. Thus, we need to know the typical loop length of the 
BP system observed bv IChandrashekhar et al.l (120131) . which can then be used as an input 
parameter in the loop simul ation. To that end, we make u se of the HMI images corresponding 
to different BPs studied by IChandrashekhar et ahl (l2013h to produce the pote ntial held lines 


of force-free magnetic helds (lAlissandrakis 


1981 


Nakagawa fc Raadd Il972l ) . A graphical 


depiction of this potential held extrapolation for a typical BP is shown in FigHJ Here, white 
lines are the potential held lines drawn using a consta nt a force-free magneti c held (with 
a = 0) connecting diherent polarities of BP7 (hg. 1 of IChandrashekhar et al.l (j2013|)). As 
these lines are very short and closes up sufficiently below the source surface, it is safe to 
use the potential held extrapolation. We calculated lengths of diherent held lines composing 
BP7. The procedure is repeated for BP2, BP3, BP5, BPS, and BPIO. BPs which had 
their opposite polarities very nearby were left out of this calculation since their loop lengths 
cannot be considered as typical. The resulting histogram of the lengths of the held lines 
is shown in Fig. [2l We observe that the histogram of lengths peaks around 10 Mm. This 
length have been used in our subsequent multithreaded nanohare heated hydrodynamic loop 
simulation. 
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Fig. 1.— Graphical depiction of potential field extrapolation for BP7. 
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Fig. 2.— Histogram of the extrapolated field lines for the BPs. 
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3. Numerical simulation of the multistranded loop 


The present simulation was per formed with a ID hydrodynamical code based on Lagrange- 
Remap method flArber et ahll200ll) . In this simulation, we constructed a global loop using 
125 strands. Each strand was evolved hydrodynamically independently of each other by 
solving the mass balance, momentum balance and energy balance equations in ID given as 
follows: 

Dp 
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Dv dp „ d'^v 

R rr. 

P = -pT . 

P 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Here t is time, s represents the coordinate along the strand which is assumed to be semi¬ 
circular, p, V, p, n and T are mass density, bulk velocity, pressure, particle density and 
temperature of the plasma respectively, and 

D _ d d 
Dt dt~^^ds' 

We require —L < s < L with L = 5 Mm, so that total length of the loop is 10 Mm. In 
equation ([2]), geos 9 represents the projected component of the solar surface gravitational 
acceleration along the semicircular loop, g being the magnitude of the acceleration. Because 
we are considering a small loop, we assume p to be a constant with a value equal to the 
surface value of 2.74 x 10^ cm s“^. Also in the same equation, the coefficient of viscosity 
u = 2.0 X 10 ^^ cm^ s ^ is assumed to be uniform throughout the plasma. In equation ([3]) , 
7 = 5/3 is the adiabatic index of the medium, while k = 9.2 x erg s“^cm“^K“^ 

is the thermal conductivity of the plasma along the loop. Furthermore, the plasma is also 
cooled th rough radiation whi ch is modeled by an optically thin radiative loss function Q{T) 
based on iRosner et ahl fjl978l) . As well, H{s,t) represents the heating function (whose char¬ 
acteristics are described in more details below) and emulates nanoflare like heating events. 
Finally in equation (H]), R = 8.3 x 10^ erg moD^K”^ is the molecular gas constant and 
p = 0.6 mol“^ is the mean molecular weight of the plasma. 

In the present simulation, t he initial and boundary conditions imposed on a strand are 
explained in equations (6)-(8) of ISarkar fc WalshI (120081) . These conditions emulate a strand 
with zero initial velocity and an initial temperature of 10^ K. The pressure and density 
initially maintain exponential profiles to represent gravitationally stratified plasma. 
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As already mentioned, the individual strands are heated due to artihcial nanoflare like 
heating events. These nanoflares can occur only in the coronal part of the loop. In the 
simulation, we chose the location of these nanoflares randomly using pseudo random number 
generator. Each such event lasted from 50s to 150s, containing energy around 10^^ ergs. The 
eff ect of a sample nanofl are containing heat 1.049 x 10^^ erg on such a strand was tested 


m 


Sarkar Sz WalshI (120081) and the result was demonstrated in Figures 1 and 2 of the same 


article. 

Due to lack of observational statistics, frequency distribution of nanoflares are still not 
well determined. If / is the frequency distribution and E is the energy, the distribution can 
be represented by df/dE = EqE~^, where Eq is a constant and fd is called the power-law 


index. In case o f flares, w 


lere s ufficient data are available the value of fd is equal to 1.8. On 


the other hand, iHudsonl fjl99ll) pointed out that to keep the corona heated predominantly 


by nanoflares the energy distribution slope has to be steeper {(d > 2). Recent observation 
of picoflares, with estimated energy ~ 10^^ erg, also suggests the power law index jd to 
be in between 2.2 to 2.7. In t he following simulations, we will primarily use fd = 2.3. 


However, I Sarkar fc WalshI (120081) has demonstrated the effect of changing fd on the global 
loop temperature, and in the present work we also verified how changing the value of fd can 
affect the overall outcome. 


The temperature profile of a loop is known to be affected depending on whether it is 
footpoint heated, looptop heated or uniformly heated (jSarkar fc WalshI 120091) . For a uni¬ 
formly heated loop nanoflares can occur anywhere along the loop. For a footpoint heated 
loop more number of nanoflares are generated close to the footpoints. Finally, for looptop 
heated loops more weight is given to the looptop. A location histogram of the above three 
cases are shown in Fig. [3l Part of our study, as described below, was dedicated to understand 
the dependence of the oscillation periods of the BP intensities on these location preferential 
heating. 


The energy input rate per unit area per unit time would also change the loop tempera¬ 
ture, thereby potentially affecting the intensity oscillation periods of the BPs. In this work, _ 

we ha ve considered two cases of fixed heating rates, one at 4 x 10® ergs cm“^ sec“^ (jSarkar fc Walsh 


20081) and another at 8 x 10® ergs cm ^ sec 


The data synthesized through the above simulation next needs to be processed to pro¬ 
duce observables (like intensities) to correspond to measurements of the AIA. We describe 
this in the following section. 
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Position in Mm 

Fig. 3.— Histogram showing preferential heating for three different cases: footpoint, looptop 
and uniformly heated loops. 




































4. Production of observables 


From the simulated data, the calculated temperature and density of a strand was folded 
through the response function of three AIA passbands, namely 171, 193 and 211 to produce 
the synthetic emission. This gives the intensity Ja,* of the strand i at passband A as follows 

h,i = Gx{T)n^{s,t)ds , (5) 


where Gx{T) is the temperature response function at passband A, Uj is the density of the 
strand which is a function of space (s) and time (t), and ds is a line element along the 
strand. One can then derive the overall intensity I\ of the global loop at passband A by 
summing over the emission from all the individual strands as follows 


125 




( 6 ) 


2=1 


In order to match with the observation of IChandrashekhar et al.l (120131) . we next needed to 
degrade the resolution of the computed intensity (l6i). For a meaningful cornparis on of the 
synthetic intensity fluctuation with that observed by IChandrashekhar et al.l (120131). we per¬ 
forme d a wavelet analysis on the aforementioned synthetic data following IChandrashekhar et al. 

(Eli). 


Like in the observation (IChandrashekhar et al.l l2013l) . we considered an eighty (80) 
minute time windows of the synthetic intensity from every simulation we have performed. 
We then convolve d the time sequence of the intensity with the Morlet function to perform the 
wavelet analysis (iTorrence fc Compel 119981) . Periodicities were calculated by removing the 
background trend by subtracting 30 minute equivalent number of points from the intensity. 
The resulting periodicities are tabulated in Table [U 


5. Analysis 


Fig. Hlshowcases the results of a wavelet analysis of a typical synthetic loop intensity. The 
top panel shows the intensity time series at the looptop. In the middle left panel we show the 
wavelet power spectrum of the same intensity time series. The cross-hatched portion of this 
panel corresponds to locations where estimates of the oscillatio n period become unreliable . 
This region is known as the cone-of-influence (COI) according to lTorrence fc Compol (119981) . 
The location of the maximum power is indicated by the white line in this plot. Next, the 
maximum measurable period is shown by a dashed line in the global wavelet spectrum plotted 
in the middle right panel. Finally, the bottom panel shows the variation of the probability 
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estimate calculated using the randomization technique associated with the maximum power 
at each time in the wavelet power spectrum. The oscillations with a probability greater 
than 95% are considered to be signihcant. The periods of the peak power are printed above 
the global wavelet plot. Here, PI refers to the most significant peak and P2 refers to the 
second most significant peak with confidence level higher than 95%. The result shows that 
the power is more for the period PI = 26.26 mins oscillation, where as the second peak is at 
P2 = 11.12 mins. 

The overall simulation results are summarized in Table [H There appears to b e no dis¬ 
tinct difference between observed oscillation periods fjChandrashekhar et al.l 120131) and the 
oscillation periods from the synthetic data of 10 Mm loops. Both the simulation and obser¬ 
vation show typical oscillation periods of around 15 - 25 mins. The characteristic oscillation 
period of the simulated loops are also seemingly independent of the heating location or the 
heat input. 

To obtain further insight into the origin of such intensity oscillations, it is instructive 
to analyze individual strand behavior. We randomly selected a strand for this purpose, e.g. 
strand no. 12 of the footpoint heated loop. Since intensity of this strand is proportional to 
the square of its density, we performed a wavelet analysis of its density squared exactly in 
the same way as have been done in the case of the global loop intensity. These results are 
also tabulated in Table [Hand we hnd similar range for the periods (15 - 25 min) as that of 
the intensity of the global loop. 


In the past. lWang et ah! (1200211 have observed oscillations of periods ranging from 14 -18 
minutes in flaring loops. It was thought that this kind of loop osci llations were excit e d by a 
single flarelike i r npuls i ve heating ev e nt. U sing ID loop simulations. iNakariakov et ah! (120041 1: 
Tsiklauri et ah! (1200411 : ISelwa et ah! (1200511 have later confirmed that indeed such oscillations 
with such periodicities can be generated using impulsive heating events. These studies have 
further concluded that these oscillations are due to slow standing waves in the loop. 

Compared to flaring loops, BP loops are signihcantly smaller in length, with the majority 
of them having length around 10 Mm as reported earlier in Section [2J The slow wave 
oscillation period P of such waves can be calculated using 


P = 


2L 

nC. 


( 7 ) 


where L is the loop length, Cg i s the sound speed o f the loop plasma, and the integer n labels 


the harmonics. Cg is given by (lAschwandenI [200411 


IMK 


[Km/s] . 


( 8 ) 























n u U-I-* ro 7 \ Noromllzed Intensity 

Probability (%) a ■ a r ■ \ 

Period (minj i | | 


10 
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Fig. 4.— Wavelet analysis of the synthetic intensity (AIA 171) of a footpoint heated loop 
when the heat inpnt rate is 4 x 10^^ ergs cm”^ sec“^. The analysis reveals PI = 26.26 mins 
and P2 = 11.12 mins. 
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Fig. 5.— Evolution of pp (in appropriate units) along the strand from the simulation of a 

uniformly heated loop with total heat 4 x 10^^ ergs cm“^ sec“^. 
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AIA 171 bp3 Light curve — Original 



AIA 171 bp3 Light curve — After trend subtraction 
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Fig. 6.— Wavelet analysis of BPS intensity (171 A) with 9 mins rnnning average snbtraction, 
revealing periods of the fnndamental mode of the standing waves. 
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Fig. 7.— Uniformly heated loop strand velocity: (top) contour plot of the strand velocity 
evolution, (bottom) velocity evolution at the center of the loop, marked by the dark line in 
the top hgure. 
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For our simulated strands with L = 10 Mm and temperature T varying between 0.2 MK 
to 2 MK, the fundamental period (n = 1), according to equations ([7]) and ([8]), should lie 
within the range of 5 to 2 mins. On the other ha nd, a direct analysis of t he syn thetic data (as 
shown in Tabled]) as well as the observations of IChandrashekhar et al.l (120131 ) yield an order 
of magnitude higher periodicity (in the range of 15 - 25 mins) for the modes corresponding to 
the two most signihcant power peaks. Clearly, the speeds of these higher periodicity modes 
are not described by equation (JHj). Indeed, e quation (IHI) provides the expression for sound 
speed of a mode which respects adiabaticity (ILandau fc Lifshitslll959l) . But the right hand 
side of the energy equation ((3|) shows a strict violation of adiabaticity over the period of the 
simulation, and one may identify the nano-flares contributing to H{s,t) being the primary 
source of this. This may be conhrmed in the plot of Fig. [5] which captures a variation of 
pp~'^ along the loop during a part of its evolution. 


Nevertheless, modes with characteristic time scales comparable with the time scale 
associated with the nanoflares should approximately respect adiabaticity. Given that the 
nanoflare lifetime varies between 50 - 150 secs, one may then identify the modes with pe¬ 
riodicity between 2-5 mins with such adiabatic modes. To uncover these modes from the 
synthetic data, an analysis was performed on the synthetic intensity (193 A) time series data 
with 9 mins background subtraction. One may contrast this with the previous wavelet anal¬ 
ysis of the intensity time series data where a 30 minutes equivalent number of background 
points were subtracted. For the present analysis, it was found PI = 3.03 mins and P2 = 2.14 
mins in case of footpoint heated loop when the heat input is 4 x 10^^ ergs cm“^ s“^. 


We next performed a similar wavelet analysis on the data of IChandrashekhar et ah 


(120131) to compare with the smaller background subtraction on the synthetic data. In par¬ 
ticular, we selected 20 mins intensity data of BP3 for this purpose. Instead of 30 mins 
background subtraction, a 9 mins background subtraction was performed during the wavelet 
transform. The result is shown in the Fig. El We hnd through this analysis that the pe¬ 
riodicities of the BP are PI = 6.1 mins and P2 = 1.2 mins, which are in the same range 
of the periodicities obtained from the analyzed synthetic data. This strongly indicates that 
the oscillations in the range of 15 - 25 mins are due to non-adiabatic modes, which either 
gets generate d by the nanoflares d isplacing the local pl asma, or through the chromospheric 
evaporation (IDoschek et al.lll980l: iFeldman et al.lll980l) due to which the plasma gets into 
the coronal part of the loop. 


Since the oscillations in the range of 15 - 25 mins, observed by IChandrashekhar et ah 


(120131) and also captured in the present simulation, are of non-adiabatic origin, one cannot use 
equation (JHj) to determine their speed. However, if one assumes that these are fundamental 
modes, then the speed of such modes turn out to be about 17 Km sec“^ according to 
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equation ([7]). Interestingly, an analysis of the general plasma flow for a typical strand (strand 
13) of the uniformly heated loop where the heat input rate is 4 x 10^^ ergs cm“^ s“^ reveals 
a velocity time evolution along the strand as shown in Fig. [71 According to this plot, the 
plasma velocity at the center of the strand can reach up to 20 Km sec“^, which is very close 
to the estimated speed of the modes presented above. Therefore, its seems highly plausible 
that the oscillations in the range of 15 - 25 mins are born as a result of hydrodynamic shock 
waves generated by the nanoflares. 


6. Summary and Discussion 


Coronal BPs have been observed since long (IVaiana et al.lll973h and are known to be 
small loop systems. It has also been known tha t the intensity of these BPs oscillate within 


certa in range of periodicities flTian et ah! [20081: [Kumar et al.l [20111: [Chandrashekhar et ah 


20131 ). However the origin of these oscillations is still unknown. In this work, we propose 


a mechanism behind such oscillation by modeling the BPs as strands of plasma heated by 
swarms of nanoflares. 


The typical loop length of such BP loops were computed using potential held extrap¬ 
olation (Fig. [1]) and was found to be aro und 10 Mm. These BP loops were then modeled 
using a ID multithreaded numerical code fISarkar fc Walsh[[2008[) . with the length set to 10 
Mm. Strands of these simulated loops were heated by series of nanohare like heating events 
during their lifetime. We considered three kinds of loops depending on the heating location, 
namely footpoint heated, looptop heated, and uniformly heated. 


The oscillation periods and temperature of the simulated loops are given in the Table [H 
We observed that diherences in the heating location do not ahect the oscillation periods. We 
also changed the heat input rate from 4 x 10^^ erg cm“^s“^ to 8 x 10^^ erg cm“^s“^ in order 
to study the ehect of increase of temperature of the loop on the periodicity. However, no 
signihcant correlation between the loop temperature and the oscillation period was observed. 
It is also worth checking if the nanoflare distribution function has any effect on the oscillation 
period. For that we ran a simulation of the uniformly heated global loop (heat input rate 
4 X 10^'^ erg cm“^s“^) with power law index (3 = 3.29. The wavelet analysis of the loop 
intensity (AIA 193) yields PI = 18.71 mins and P2 = 9.35 mins, when 30 mins background 
substruction was performed. This suggests the observed oscillation period is unaltered with 
the change of power law index (3. 


Table [T] shows that the oscillation periods of the simulated loops were in the range of 15 


25 mins and such range compares well with the values observed by [Chandrashekhar et ah 
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(120131) ■ These oscillations are present both in the intensity time series of the global loop 
as well as in the flnctnations of the individnal strand density. It is snggested that these 
oscillations are dne to non-adiabatic shock waves present in snch mnltiple nanoflare heated 
loops, since their periods are mnch higher than those for adiabatic modes generated for the 
temperatnre range 0.2 - 2.0 MK of the loops, and also becanse the speed of these modes are 
in the same range as the typical speed of the plasma (~ 20 Km sec“^). In fact, for the given 
temperatnre range of the plasma, the adiabatic modes have periodicities in the range 2 - 
5 mins according to eqnation ([H]). These adiaba tic modes can be nncovered b oth from the 
simnlated data as well as from the observations fjChandrashekhar et al.l (120131) . from the in¬ 


tensity light cnrve of BP3, a typical BP) throngh a 3 mins bac kgronnd snbtraction in s tead o f 
a 30 mins backgronnd snbtraction as originally performed by Chandrashekhar et ah (j2013 ). 


To snmmarize, in this work we have proposed two kinds of oscillations in coronal BP 
intensities, namely oscillations dne to shock waves and oscillations dne to adiabatic modes. 
Farther conhrmation of this proposal can be achieved by stndying the Doppler velocity map 
or non-thermal line broadening of BPs. These velocity maps and line broadenings can then 
be compared with the simulated synthetic data resulting in a better understanding of such 
waves. 


We are grateful to the anonymous referee for carefully reading the manuscript and 
suggesting us important changes. The Center of Excellence in Space Science India (CESSI) 
is supported by the Ministry of Human Resource Development, Government of India. 








Table 1: Periodicities observed in simulated loops. 



Energy input 



Intensity 



Density squared 

Looptop temperature 



171 

193 

211 

(single strand) 

(MK) 


{ergs / cm? / sec) 

PI 

(mins) 

P2 

(mins) 

PI 

(mins) 

P2 

(mins) 

PI 

(mins) 

P2 

(mins) 

PI 

(mins) 

P2 

(mins) 


Looptop 

4 X 10^7 

15.7 

24.3 

17.15 

10.2 

17.2 

9.35 

22.3 

8.6 

2.5 

heated 

8 X 10^^ 

- 

- 

22.3 

11.1 

20.4 

11.1 

25.3 

11.1 

3.1 

Footpoint 

4 X 10^7 

26.3 

11.1 

11.1 

2.6 

11.1 

26.5 

17.1 

13.2 

2.2 

heated 

8 X 10^7 

- 

- 

22.3 

7.2 

20.4 

7.2 

14.4 

7.9 

2.7 

Uniformly 

heated 

4 X 10^7 

17.2 

26.5 

17.2 

7.9 

- 

- 

18.7 

8.6 

2.2 
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